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Artificial Intelligence Device and Corresponding Methods for Selecting 

Machinabilitv Data 

Field of Invention 

5 

The invention relates to an artificial intelligence device and corresponding methods 
for selecting optimal machining parameters especially for cutting speed and feed 
rate control, and more particularly applicable to computer controlled milling 
machines, drilling machines, grinding machines, turning machining and other such 
10 machines. 

Background of Invention 

Machinability data is never defined preciously in a scientific way. Machinability 
15 data consists of the selection of the appropriate cutting tools and machining 
parameters, which i ncludes cutting speed, feed rate a nd depth of cut. It plays an 
important role in the efficient utilization of machine tools and thus significantly 
influences the overall manufacturing costs. It has become more and more important 
especially in computer aided design and computer aided manufacturing. 

20 

The main objective of machining is to satisfy the demand of the work piece 
dimensions and its surface finish. The machinist is subjected to a wide range of 
possible speed, feed rate, depth of cut, cutting tool selection, material work piece 
and other machining parameters. The application of machining economics to shop 
25 floor practice has been 1 imited a nd the correct choice of machining conditions is 
normally overlooked. 

The general practice o f e xperienced m achinists i s captured, f ormalized and made 
available in the form of machining data handbook, data sheets and charts, 
30 nomograms and slide rules. Metcut's Machining Data Handbook (MDH) is among 
the common sources. The data in the handbook were separated into different types 
of machining process, such as face milling, end milling, turning, boring, etc. The 
data were grouped according to the types of work piece material and the respective 
hardness. The handbook provides the speed of cut and the feed rate corresponding to 



2 



the depth of cut. For each work piece material and hardness, one may scan through 
different possible tool materials. Then, speed and feed rate are selected according to 
tool-work-piece material combination, depth of cut, and finishing condition. 
Although handbook approach is often a logical and effective solution to the 
5 requirements of machinability data, it has the following limitations: 

• Handbook recommendations represent a starting set of cutting conditions and 
hence tend to be conservative in order to cope with a worst-case machining 
scenario. 

• Handbook data only applies to a particular machining situation. This data 
10 may not be suitable for slightly different machining situation. 

• Handbooks are manually input-output oriented, and hence lack compatibility 
with the objective of integrated automation of the manufacturing system. 



Artificial intelligence (AI) is normally referred as any operation being carried out 
15 based on the wistful hope that at least some of the flexibility and power of the 
human brain can be reproduced by artificial means. The characteristics of brain 
function have inspired the development of artificial neural network. The way of 
human making decision and analyzing in real life has inspired the development of 
fuzzy logic. The experience of the machining operators may express in terms of 
20 fuzzy logic and/or modeled with artificial neural networks. 



Turning operation of wrought carbon steel with different cutting tools has been 
modeled successfully with artificial neural networks and the fuzzy reasoning have 
shown promising results in prior-arts by the inventor. Metcut's Machining Data 
25 Handbook was chosen as the reference for validation and data source training. The 
methods are further generalized in the present invention to cover all machining 
process with different work piece materials and cutting tools. 



Summary of Invention 

30 

Accordingly, it is a first object of the present invention to provide an artificial 
intelligence device and corresponding artificial intelligence methods which provides 
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a better solution compared to traditional methods for selecting optimal machinability 
data. 

The present invention comprises of the three main embodiments. The proposed 
5 systems works with different or combination of different AI methods. Among them 
are Fuzzy Logic, Neural Networks and Genetic Algorithm. The embodiment will 
carry different function for different AI models. 

In a first aspect of the present invention, the device comprises of Input Component, 
Inference Engine, and Output Component. The Input Component comprises of 

10 different input interfaces, means to store and manipulate the collected inputs. The 
corresponding input manipulation depends on the method incorporated in the 
inference engine. The Inference Engine will process the input and prescribe the 
outputs in respective form depending on the AI model incorporated. The yielded 
outputs, which are the recommended initial cutting parameters, will be further 

1 5 processed by the Output Component. The output will then be further processed to 
turn into required form and format. The transformed output will then be fed into the 
machine for further action by the machine. 

Brief Description of the Drawings 

20 

Figure 1 is a block diagram showing the Machinability Data Selection System with 
Artificial Intelligence. 

Figure 2 is the inputs fuzzy expression. 

25 

Figure 3 is the output fuzzy expression for cutting speed fuzzy model. 

Figure 4 is the min-max range of fuzzy membership functions for cutting speed 
fuzzy model. 

30 

Figure 5 is the Fuzzy rules for high-speed steel fuzzy model. 
Figure 6 is the Fuzzy rules for uncoated brazed carbide fuzzy model. 
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Figure 7 is the Fuzzy rules for uncoated indexable carbide fuzzy model. 

Figure 8 is the Fuzzy rules for coated carbide fuzzy model. 

5 Figure 9 is the input and output membership functions for high-speed steel feed rate 
fuzzy model. 

Figure 10 is a flow chart for general system of fuzzy rules genetic optimization. 

10 Figure 1 1 is a block diagram of simple non-linear neural network of machinability 
data for turning 

Detailed Description of the Invention 

15 With reference to Figure 1, the basic aspect of this invention is the use of Artificial 
Intelligence techniques for selecting optimal machining parameters in conventional 
machining operations. The device comprising three main components which are the 
input, the inference engine and the output. 

20 The inputs are the machining operations, work piece material, machining tool type 
and depth of cut. The outputs are the machining parameters comprising of the 
optimal cutting speed and feed rate. The inference engine can be established using 
AI techniques such as fuzzy logic, neural network or fuzzy-neural network. 

25 In the first embodiment, there are three basic components for a typical fuzzy expert 
system which are input fuzzification, rules application (or inference engine) and 
output defuzzification. The input fuzzification translates the system-input variables 
into universe of input memberships. There are two common methods used for an 
inference engine which are Max-Product method and Max-Min inference method. 

30 The difference is the aggregation of the rules. They use multiplication and truncation 
of the output fuzzy set with the yielded result respectively. The defuzzification 
process is defined as the conversion of a fuzzy quantity, represented by a 
membership function, to a crisp or precise quantity. There are two common methods 
for defuzzification of fuzzy quantities which are the Max method and Centroid 
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method. Two separate models have been developed to yield the desired output for 
wrought carbon steel turning process which are the Cutting Speed Fuzzy Model and 
the Feed Rate Fuzzy Model. 

5 For the Cutting Speed Fuzzy Model, the hardness of work piece material and the 
depth of cut are the inputs and cutting speed is the output of the fuzzy model. Figure 
2 and 3 show the fuzzy expressions for work piece material hardness, depth of cut 
and cutting speed respectively. Figure 4 shows the range of the fuzzy input 
membership functions and the output membership functions. All membership 

10 functions (both inputs and outputs for all models) are of triangle shape and are well 
distributed. Four fuzzy models with different fuzzy rules amongst themselves are 
developed. The first model with fuzzy rules shown in Figure 6 was designed for 
high-speed steel tool. Other models were then designed to get optimum result 
(cutting speed) for uncoated i ndexable carbide tool, uncoated b razed carbide tool 

15 and a coated carbide tool. The rules set for these tools are shown in Figure 5 to 7. 
The search for the optimum model for every cutting tool has been limited to the 
variance in rules only, where the input and output fuzzy sets are maintained and 
linearly proportionate among two types. 

20 The Feed Rate Fuzzy Model is a single-input-single-single-output fuzzy model. The 
fuzzy shape of the input and output fuzzy sets is that of a triangle, but not of an 
equal-sided triangle. Figure 8 shows the fuzzy expression for both the input and 
output membership functions. The model applies straightforward fuzzy rules, which 
are 

25 i. If the depth of cut is very shallow then the feed rate is very slow 

ii. If the depth of cut is shallow, then the feed rate is slow 

iii. If the depth of cut is medium, then the feed rate is medium. 

iv. If the depth of cut is deep, then the feed rate is fast 

v. If the depth of cut is very deep then the feed rate is very fast 

30 Unlike the Cutting Speed Fuzzy Model, the Feed Rate Fuzzy model maintain the 
fuzzy rules as a common rule while changing the input and output membership 
functions in order to yield the best results. The high-speed steel fuzzy model is set as 
the reference membership function throughout the new membership function 
development for other tools. 
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The following equation El is used to define the input and output fuzzy sets of other 
tool type models. 



5 El 
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The corresponding new value of a membership function (input and output as well), 
Vnew is yielded from the value of the input and output membership functions of the 

10 High-Speed Fuzzy Model, V st d (so named reference fuzzy model), 

> V st d,max and V s td,min are the maximum and minimum value of the new membership 
functions and the standard membership functions. The exponential value, n, will 
determine the pattern of the yielded membership function, n = 1 will lead to a linear 
proportional relationship. To obtain a better correlation, the value of n is obtained 

15 from Equation E2 instead of being a constant. 
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The value n will be equal to A at the lowest V st d valued and gradually decrease to B 
(which must be a value smaller than unity) when V st£ i is at the required ratio of the 
25 range, R. The changing pattern of n is dependant on the constant Pi. Then n will 
gradually increase to unity at the maximum value of V stc j. Again, the changing 
pattern of n is depending on constant P 2 . 

In the second embodiment the genetic optimization is used to further optimize the 
30 fuzzy rules. Major operations of genetic algorithm depend on random choice. A 
random number generator class is developed. It consists of three main functions. 
They are: to generate a random real number from 0 to 1, to generate a random 
number from a user given start to a user given end integer value and to give green or 
red light for a user given probability. 
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Genetic optimization of fuzzy rules has been carried out with the help of an object 
oriented genetic optimization library (GOL). It consists of several useful and inter- 
connected classes. To evaluate the fuzzy result, the Fuzzy Set Handling (FSH) class 
5 has to be included. Figure 10 shows a general system flow of the genetic 
optimization process for fuzzy rules design. 

The population needs to be initialized. Information such as the number of alleles and 
the length of each allele are predetermined. The GOL use bit-wise interpretation, 

10 which means, the length of a particular allele is expressed in term of bit. If an allele 
carries a possible value from 0 to 7 (or 8 possible features), the length of allele is 3. 
The fuzzy models consist of 5 fuzzy sets for all inputs and 15 fuzzy output sets. The 
system required 25 fuzzy rules with 15 possibilities each. For initialization, 25 
alleles are required and the length of each allele will be 4 bits to cope with 15 

15 possible values. Besides initialization of the alleles, the probabilities of crossover 
and mutation, and the size of population are needed. Proper size of the population 
and probabilities will yield better optimization results in terms of speed. 
The population consists of 80 sets of fuzzy rules. The crossover and mutation 
probabilities are set as 0.6 and 0.009, respectively. The initial members of the 

20 population are required. The fuzzy rules from Fuzzy Logic are assigned as one of 
the initial members. The rest of the rules are generated automatically and randomly. 
Initialization of FSH class in Figure 1, referring to the construction of the FSH class 
includes building then membership functions of inputs and output. With the FSH 
class, the fitness of each member of the initial population has to be calculated. The 

25 one with better fitness stands a better chance of success in reproduction. In the 
process, crossover and mutations may occur and thus, create new species. Crossover 
operation is single-point operation. The mutation operation is a random change in an 
allele, instead of a bit in the allele. The new generation has a tendency towards 
better quality. The reproduction continues until the population is completed. The 

30 whole process repeats until the maximum number of generation is obtained. Fitness 
is obtained through Equation E3. 



E3 
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The value of j in E3 is the numbering of the predetermined inputs starts from 0. 
Thus, value n is the ending number of the predetermined inputs, in this case 79. 
Error factor is for rules violation penalty. The value K is arbitrarily selected positive 
5 value. It must be foreseen to be more than then achievable maximum summation of 
the absolute error percentage mean and the error factors. Generally, with lower mean 
absolute error percentage, the fitness value is higher and nearer to K. For the final 
consideration of the fitness in the competition of reproduction selection, Equation 
E4 is employed. 

10 

E4 FFitness, = Fitness i — m'm(Fitness 0 ...Fitness n ) 

The Genetic Optimization Class allows control of fuzzy rules pattern, which means 
that the designer can provide his/her expertise in specifying the relationship between 

15 the inputs and the outputs. Each violation will cause a penalty value, thus more 
violations will cause lower fitness. The user can set the penalty value. The penalty 
value must not be too small, which will not cause a significant effect. It must not be 
too large either, which will make the whole optimization process inefficient. The 
summation of all penalty values must not be greater than the difference between 

20 value of K and the mean of all absolute error percentage. 

In the third embodiment, a non-linear artificial neural network, which consists of 
both the novel product neurons developed by the inventor and the common 
summation neurons, is designed to prescribe machinability data for turning process 

25 of wrought carbon steel with different tools. Two feed-forward neural networks are 
suggested. Figure 11 shows the artificial neural network suggested. The network 
incorporates simple non-linearity of second input (the depth of cut). Two different 
neurons are indicated with different symbols, with "E" representing normal 
summation neuron and "IT" representing the new Product Neuron. Transfer function 

30 of respective neuron is placed following the neuron symbol. "PL" stands for pure 
linear transfer function and C P2' stands for squared transfer function. The design of 
the network depends on the relationships among the parameters. The artificial neural 



networks suggested are supervised feed forward networks. No closed synaptic loop 
is used. 

Summation neuron is the classical neuron, which is used in most of the present 
5 neural systems. A neuron with R inputs and bias, b has the network summation 
value, n computed by the following equation: 

n = Wp + b 

where W and p are the weight matrix and input matrix, which store corresponding 
10 weights and input values respectively. The transfer function, f may be a linear or 
non-linear function of n. The output neuron, a, can be written as below: 

a=/n) 
=/Wp+b) 

15 Product neuron is introduced by the inventors to deal with non-linearity and 
establish cross-neuron relations in a simpler way. The neuron may be hooked with 
multiple neurons but only weighted once. Similar to summation neuron, it holds a 
bias as an independent variable. The n is calculated with Equation E5. 



20 E5 n 
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To avoid complication in design of learning methodology and data storage, the 
neuron is deliberately planned to have one overall weight, instead of individual 
weight of each input. Thus, processing (including learning) of a single Product 
25 Neuron is simpler than processing a summation neuron. 

The final output(s) of the network are computed with propagation forward through 
the network from the input(s). The forward propagation process can be expressed 
with the Equations 6, 7 and 8, where a m is the output matrix of all neurons in layer 
30 m. The initial matrix a 0 represents the input values. The outputs of the whole 
network, which has M layers, are gathered as matrix a. The matrix f" 1 collects 
transfer functions of all neurons in layer w. 
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The back propagation algorithm is a generalization of the least mean squared 
method to provide learning ability to the network. It is capable of handling multi 
layer network. The classical back propagation algorithm is summarized in this 
section. 

The first step of the algorithm is to find the final output of the network with forward 
propagation. Then the sensitivity of each neuron is calculated with backward 
propagation through the network from the outermost layer, which yields the 
output(s). Recurrence relation is established among the neurons and chain rules are 
used to obtain the sensitivities. The neuron sensitivities of last layer are calculated 
Equation E9 in matrix form. 
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Matrix t consists of all target outputs. The F M (n M ) is the sensitiviy function of 
respective neuron and is expressed as the Equation 10, where Sm is the number of 
neurons in layer m. For other layers, the Equation 1 1 is used. 
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Finally the weights and biases are updated with selected update method, which may 
5 be referred as learning algorithm. 

Weights and biases are updated with Approximate Steepest Descent method in this 
paper. The weights and biases are updated as shown in Equations 12 and 13. 

10 E12 W m (k + \) = W m (k)-as m (a m - l y 

E13 b m (k + \) = b m {k)-as m 

Momentum is used to smooth out the oscillations in the trajectory towards the 
15 optimum location. Equations 14 and 15 are used to achieve such task, where AW m 
(k) and Ab m (k) are the Ath change of weights and biases at layer m. The learning rate 
and the momentum coefficient are denoted as a and y respectively. 
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E 1 4 AJV m (k) = yAW m {k - 1) - (1 - y)as m (a m ' l ) T 

El 5 Ab m (k) = yAb m (£-!)-(!- y)as m 



From the mathematical point of view, the BP learning with SD method only 
promises a general minima achievement, may it be a global, local or weak minima. 
25 It is guaranteed to converge to a solution that minimizes the mean squared error, as 
long as the learning rate is not too large. Typically, the initial weights and biases are 
chosen to be small random values. In this way, one can stay away from a possible 
saddle point at the origin. It is also useful to try several different guesses in order to 
be sure that the algorithm converges to a global minimum point. 



A new retraining scheme is suggested to move the solution away from any local or 
weak minima during the training of an artificial neural . network. If the network is 
found having weak response towards a particular training set, extra retraining are 
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enforced until the network's performance index of the training set is under a 
predefined index of the training set is under a predefined-level, a. Generally, 
implementation of the scheme incurs certain level of interruption towards the 
network (refer to the weight and bias values). Thus, one may still observe a sudden 
5 fluctuation on the performance surface even the solution reached its global minima. 
This can be solved by incorporating a simple global-minima tracking procedure. 
Implementation of the RR does not require Variable Learning Rate (VLR) 
algorithm. The involvement of the VLR will slow the overall efficiency of the RR 
scheme, as the main purpose of the RR scheme is to bring the solution out of the loal 

10 or weak minima towards the global minima. Initial studies by the inventors show 
that the VLR with RR cannot bring the solution out of the local or weak minima 
efficiently. The RR scheme requires a stable learning rate to obtain overall effect 
during the training process. The SD method alone is good enough to carry the task. 
The squared error limit, a, is calculated with Equation 16, where ajnitiai and a en d are 

15 the user defined values which indicate the squared error limit at the initial and the 
end of the training respectively. iTrain and iTrain max are the current number of 
training and the maximum number of training desired respectively. As general 
practice, the value of ajnitiai should be larger than a en d- 

20 E16 g = ainitial ~ ^ end - x iTrain 

iTrain max 



The scope of the invention should be defined only in accordance with the claims that 
follow. In the following claims, reference characters used to designate claim steps 
25 are provided for convenience of description only, and are not intended to imply any 
particular order for performing the steps. 



